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The 4 sites and 8 sites ID anti-ferromagnetic Heisenberg chains in the Jordan-Wigner representa- 
tion are investigated within the standard Hartree-Fock and RPA approaches, both in the symmetry 
unbroken and in the symmetry broken phases. A translation invariant groundstate, obtained by the 
projection method as a linear combination of a symmetry-broken HF state and its image under re- 
flection, is also considered, for each chain type. It is found that the projection method considerably 
improves the HF treatment for instance as far as the groundstate energy is concerned, but also with 
respect to the RPA energies. The results are furthermore confronted with the ones obtained within 
so-called SCRPA scheme. 

PACS numbers: 75.40.Gb, 75.50.Ee, 75.10.Pq, 75.10.Jm 



I. INTRODUCTION 

The ID Anti-Ferromagnetic Heisenberg Model (AFHM) Q is a prominent many body research field. It was the 
first many-body model to be solved exactly by Bethe with his famous ansatz || and, as a solvable model, it has 
proved to be very useful for testing and developing many body approaches ever since. In spite of tremendous progress 
in the understanding of many aspects of the model, it remains a very active field of interest and research [HQ- 

In [f| we have considered the anisotropic Heisenberg model 

N f 1 + + z\ 

in the Jordan-Wigner representation. We have found that for low values of the anisotropy factor g, the Hartree-Fock 
(HF) groundstate is invariant under translations. However, for large g values a (spherical) HF state is no longer 
stable. Instead one has to introduce a symmetry broken (deformed) HF state which becomes degenerate. In a 
non translational- invariant HF groundstate was introduced in order to avoid unstable RPA modes, i. e., RPA modes 
with purely imaginary or complex energies. However, it is well known that the exact groundstate of the Heisenberg 
model is not degenerate (except for g = oo), so that the occurrence of a degenerate HF groundstate is a spurious by 
product of the HF approximation. It is therefore natural to look for an improved description by considering linear 
combinations of deformed HF states. This is the main objective of the present note. We are therefore lead to introduce 
a projection method, implementing the Peierls-Yoccoz projection approach in the version variation after projection 
Q. As a result, considerable improvement of the HF treatment is achieved, not only as far as the groundstate energy 
is concerned, but also with respect to the RPA energies. We also give the results obtained with an extension of RPA, 
the so-called SCRPA (self-consistent RPA)@. It may also be observed that the projection method results become 
exact for the 4-spin AFHM. 
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II. FORMALISM 



A. Variation after projection method 

We consider the one dimensional Heisenberg antiferromagnct model pj in the Jordan- Wigner representation [?} ■ In 
momentum space, the Hamiltonian may be written [f| 

H =^+J2 £ °q^1 + ^ E ^1529354^4 V>94 %3 
q 9192 93 94 

where the ipt > "09 are the usual fermion creation and destruction operators with momentum q and 



e q = cosq - 1, V qiq2q3qi = — (cos(gi - q 3 ) + cos(g 2 - qi))S(qi + q 2 - 93 - <&) 



N 
with 

_ Jtt f j = ±1, ±3, ■ • • ± (JV — 1), ifiV/2iseven 

9 ~ AT ' \j = 0, ±2, ±4, • • • ± (JV — 2), N, if iV/2 is odd 

We concentrate on chains with equal numbers of spins up and down, so that N/2 is the number of fermions. The 
Hartree-Fock (HF) state is obtained by filling up the lowest N/2 single particle levels. The set of these levels may be 
written 

r 7T 7T ~| 

D = jg : - < q < n or - n < q < --J , 
so that the (translation invariant) HF state reads 

i*> = n ^i°>- 

h&D 

It has been shown in Q that some frequencies of the RPA modes associated with |$), with momentum q close 
to 7r, (i.e., q > q c « 0.9l7r), are complex, implying that this state is unstable and other HF states exist which are 
energetically preferred. We consider the canonical transformation which couples levels with momenta q and q — it, 

(I = <*h1>1 + Ml-«, CL, - « h ^_ w - Ml he[^]u [-vr, -f ] = A 

where ah, Phi are rea l parameters such that + = 1. This transformation breaks the translational symmetry of 
the model. We say that the basis of the states ip q is spherical and the basis of the states ( q is deformed. It is well 
known that, very often, an energy gain may be obtained in mean field approaches by breaking some symmetry of the 
Hamiltonian. The HF state, 



m = n cli°>. ( 2 ) 



heD 



destroys translational invariance and leads to a lower energy expectation value for some non-vanishing values of the 
parameters otk, /3fe, which should be determined variationally. Indeed, let us consider the unitary translation operator 



u n = c mP , p = J2q4A- 

We have 
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Thus, |$) is an eigenstate of U n , but |$) is not. The deformed HF states |$ n ) = U n \$>), are not all distinct. There 
are only two distinct ones, for even and for odd n. Indeed, |$o) = ^2) = ■ ■ ■ = |$) and = |$3) = • • ■ = J7i|$). As 
Figs. Q] and [2] show, the groundstate HF energy, based on the state |i>), which breaks translational symmetry, greatly 
improves the prediction of the symmetrical state |$). However, in the spirit of the Peierls-Yoccoz method [6[, it is 
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Figure 1: (Color online) Four site system. Groundstate energy versus g. The green dotted line shows the non-deformed HF 
result. The blue dot-dashed line shows the deformed HF result. The red dashed lines shows the projection method result 
(variation after projection), superimposed on the exact result shown by the black solid line. 

natural to look for a still improved description in terms of a linear combination of the degenerate deformed states, 
|$) and |$i), which amounts to considering the symmetrical projection of |$). Since U\{1 + Ux)\&) = (1 + J7i)|$), 
the state (1 + E/i)|$) is translation invariant. It is an eigenstatc of P with zero eigenvalue. We describe briefly the 
computation of the expectation value of an arbitrary many-body operator X in the projected state, 

($|1 + J7i|5) 

While |$) is given by ([2]), the reflected state C/i|$) reads 

= n ^ + io) = n tficMio>. (3) 

heD heD 

We have 

UiClu} = Ch = UhTpl - M f h _ n = (a 2 h - pl)d - 2a h p h (L„, 
UrCt.M = C'U = a h ^ h _ n + Ml = (a 2 h - ft)Cl_„ + 2a h Ml 

This implies that the particle-hole concept breaks down in connection with the projected groundstate. The determi- 
nation of ($|J7i|$) and ($|Xf7i|$), where X is some operator, involves a large number of contractions with respect 
to the vacuum |0). Obviously 

wii*> = n (oiaciio) = n k - pi). 

heD heD 
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Figure 2: (Color online) The same as Fig. [T] but for the eight site system (N=8). 
The computation of ($|Xf7i|$) is also straightforward. We exemplify it for the case of a one-body operator, X = 

^ ] 2 ; pip 2 Cp 1 Cp2) 



($|XJ7i|$) = ^2 {x hh a\ + (x h{h _n) + X( h _ v ) h )a h p h + .T(, l _ 7r )(/ l _ 7r) ^, 2 l ) J| (a\ - Pi) 

heD \ h^keD 

The expectation value of the Hamiltonian in the projected state may be written 

{$\H(1 + Ui)\§) 



($|1 + C/i|$) 



(4) 



where we have explicitly indicated the dependence on the parameters ■ ■ ■ Pn/2i with respect to which minimization 
should be implemented in the end. 

In Fig. [TJ the performance of the three approaches, HF, deformed HF and variation after projection are compared, 
for the 4 spin system, in connection with the estimation of the groundstate energy. It is remarkable that the result 
obtained by the variation after projection method, coincides with the exact result. In Fig. [21 the corresponding results 
for the 8 spin system are presented. Again, the performance of the projection state method is remarkable. 



B. Random phase approximation method 



It is well known that the stability of the HF state |<f>) may be investigated by studying the perturbed state 
|$s) = exp(iS')| < I > ) where S is an infinitesimal one-body Hermitian operator. Expanding the expectation value of the 
Hamiltonian, we have 

($s\h\$ s ) = (mm + l<ns, w,s\m + • • • , 

where, in the right hand side, a term — z($|[5, has been omitted, since, by the HF condition, it vanishes. 

Stability of |$) means that (^[S 1 , [H, S]} |$) is positive semi-definite. 
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The mean field dynamics, which is prescribed by the time evolution of is determined by the quantal Lagrangian 

C = - (*s|*s» - |$s), 

where the dot sign denotes time derivative, as usual. For small S, we may replace C by its quadratic part, 

cw = l -(<s>\[s,s\*)-\(ms, [n,sm, 

the time dependence being encapsulated in S. The time evolution of S is determined by the equation 

i($|[S,5]$)-($|[H, [H,Sm = 0, 

where S denotes an arbitrary one-body operator (for instance, 2 = ip^tpg). The ansatz S = cxp(—iujt)<d leads to the 
equations of the so-called random phase approximation (RPA), 

w^p.eil^-^HH, [H,Q\m = 0. 

If to is not a real number, then |$) is not stable and should be replaced by a stable solution. This is what is done when 
the symmetrical state |$) is replaced by the deformed state |<£>) for g > g c , the critical strength for the emergence 
of the instability. Although, in the conventional RPA, |$) is a Slater determinant, extended versions of the RPA, 
for which |$) is replaced by a state containing correlations, have been proposed and may be found in the literature 
@, IE El- ln the case of the Heisenberg model in the Jordan- Wigner representation, it is advantageous to replace |4>) 
by the projected state (1 + t/i)|$), in the RPA equations. This is permissible because, for p ^ q, we have 

($|(i + tr 1 )[^, ) fl](i + i7 1 )|*) = o, 

since the projected state has vanishing momentum and the operator iptipq carries a finite momentum p — q. The 
ket |$) entering in the projected state should be fixed by the variation after projection procedure [f|, that is, the 
parameters • • • , (3n/2 should be determined by minimizing £q(Pi, • • • , Pn/2) given by 

In order to determine the RPA frequencies, it is enough to consider the ground state expectation values of the 
operators 

N {q 2 P2]Piqi) = fyt^qz^l^qi] = (V^VvWi - V^V^Wi) ( 5 ) 

T (q2P2]Piqi) = [^l^i^Sa^a^l^qA] = (e Pl ~ S qi ) (V>| 2 V'gi <W ~ V'p 1 V ; g2<W) (6) 

a 

Ov(q2P2-,Piqi) = bPl^Ao 5Z Vabcdlpllpt^dlpc^lpqA} 

abed 

= 2 Y1 ^bcp^l^t^q^c + 2^ Vq^cv^l^li'q^c 
be be 

- ^2 Va.bvivii'li'li'q^qi - V ni2od l >Pi 1 tpl 2 tpd'>Pc 

ab cd 

- 51 V qibcd1pl 2 lpilpdi'cSq 2 p 1 ~ ^ Vabpid^a^l^d^qJqiPi- (7) 
bed abd 

It is convenient to consider in above operators pi — qi = 92 — P2 what implies momentum conservation. The 
commutators and double-commutators have been worked out in the symmetrical basis, although the deformed basis 
is equally admissible for this purpose. However, the expectation values defining the RPA matrices should be taken in 
the appropriate state: the spherical HF state, the deformed HF state or the projected state, sticking to the variation 
after projection prescription. The RPA frequencies are the eigenvalues u> of the eigenvalue problem 

uN( = (T + gV)(, 

where £ denotes the eigenvector and A/", T, V are appropriate matrices. Explicitly we have, 

U •^ r (<?2P2),(pi<?l)Cpi<Ji = (^(92P2),(Pi9i) 9^(q 2 p 2 ),(piqi)) < ~>Piqi J 

Pi 91 Pl9l 
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Figure 3: (Color online) Four site system. RPA frequencies vs g. The modes for momentum transfer q = — are degenerate with 

the upper mode for the momentum transfer q — n. The RPA results are represented by the red dot-dashed lines. The results 
for the projected state and SCRPA shown by the blue dashed lines and the green dot-dashed lines, respectively, coincide with 
the exact results which are shown by the black full curves. 



where 



■ / ^(92P2),(pii3l) 

^ / '(92P2),(Pl9l) 



(<f>\0 N (q 2 p 2 ;piqi)\G>), 
($|0t(«m>2;pi«i)|S}, 

(®\O v (q 2 p 2 ;piqi)\$). 



(8) 
(9) 
(10) 



The RPA frequencies based on the projected state are the eigenvalues u of the eigenvalue problem 

u(M + M'K = ((T + V) + g(V + V'))C, 

or, explicitly, 

+ - /V "')(92P2),(pi9i)Cpigi = ^2((T + T')( q2P2 ),( piqi ) +g(V + V')(g 2 p 2 ) 1 (p l9 i))Cp I «i) 



where the matrices Af' , T', V are given by 



A/"' 

(92P2),(Pl9l) 

T' 

'(92P2),(Pl9l) 

V' 

(?2P2),(Pigi) 



(^\0 N (q 2 p 2 ;p 1 q 1 )Ui\^), 
($\0 T (q2P2;Piqi)U 1 \$>), 
{®\O v (q2P2-,Piqi)Ux\$). 



(11) 
(12) 

(13) 



The computation of the matrices Af' , T', V is straightforward, although somewhat more labor consuming than the 
computation of the matrices Af,T,V. Some computational details can be found in Appendix A. 
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Figure 4: (Color online) Eight site system. RPA frequencies vs g for the momentum transfer q = 
shown by the black full curve. The performances of S-RPA, D-RPA, P-RPA and SCRPA are compared. 



The exact results are 



Besides the evaluation of RPA with the symmetry projected state, we also show results corresponding to the so- 
called Self-Consistent RPA (SCRPA) which was applied to the AFHM in our previous publication [f|. The principle 
of this approach is explained in that reference but for completeness, let us just outline the main aspects of SCRPA. 
It is based on a particle-hole (ph) excitation operator of the form 



ph 



where the excited state is given by QjjO) = The ground state should be the vacuum to the corresponding 

destruction operator, that is Q^lO) = 0. Defining an average excitation energy f2 M via the energy weighted sum rule, 
i.e. 



a. 



[H,Qi 



]|0) 



(o\[QM]\o) 



and minimizing with respect to the amplitudes X, Y, yields RPA type of equations (0\[SQ^ 



[H,Q* 



r2 M (0|[<$Q M , Qjj|0) with the correlated ground state defined via the killing condition of above. 



Sfi) L"i "tjuj 

In the double com- 
mutator appear at most two body densities and also one body densities which can be expressed as a functional of the 
X, Y amplitudes. Self consistent equations for the amplitudes X, Y are then obtained which can be solved numerically 
by iteration. More details can be found in ref @. 




Figure 5: (Color online) Eight site system. RPA frequencies vs g for the momentum transfer q = — . The exact results are 
shown by the full curve in black color. The performances of D-RPA, P-RPA and SC-RPA are compared. 



III. EXCITATION ENERGIES: RESULTS AND DISCUSSIONS 



A. The four sites system 



Results for the four sites chain are summarized in Fig. El where we compare the excitation energies according to 
the spherical RPA (S-RPA), the deformed RPA (D-RPA) and the RPA based on the projected state (P-RPA) with 
the exact excitation spectrum. For g < 1.4, S-RPA shows some contact with the lower exact branch, but this contact 
is completely lost by D-RPA, for g > 1.4. However, the excellent performance of P-RPA, which precisely reproduces 
the exact spectrum, is remarkable. 



B. The eight sites system 



We consider, now, the case of the eight sites chain. In figures El El El El and \E\ we show the excitation spectrum as 

7T 7T 3tT 

function of the coupling constant g (anisotropy factor) , corresponding to the momentum transfer q = — , — , —— , and 

7r, comparing the results of S-RPA, D-RPA and P-RPA, with the exact ones. For g ~ 0.7, S-RPA becomes unstable. 
Nevertheless, we represent it in Figs. HI El El El f° r comparison and because in some cases it appears, unexpectedly, 
to describe better the physical situation than D-RPA. In Fig. 01 we notice that the performances of SCRPA, P-RPA 
and D-RPA, which for g < 0.7 goes over to the S-RPA, are of similar quality, although P-RPA slightly improves over 
D-RPA. However, beyond g ~ 1, D-RPA and also P-RPA depart significantly from the exact solution. Only SCRPA 
stays reasonably close for large values of g. The latter is somewhat surprising, since SCRPA is here evaluated only in 
the 'spherical' single particle basis which, in the HF approximation, stops to be valid beyond g ~ 0.7 (as will be seen 
later in Figs. [Hand El). 

One may conclude that in what concerns the momentum transfer q = ir/4, all approaches, S-RPA, D-RPA, P-RPA 
perform reasonably well below g < 1. The more reasonable agreement of SCRPA with the exact values for g > 1 may 
be accidental. Surprisingly P-RPA does not improve much on D-RPA. We will see that this is not always the case for 
the other transfers and that quantum fluctuations apparently arc of different importance in the different channels. In 




Figure 6: (Color online) Eight site system. RPA frequencies vs g for the momentum transfer q = — . The exact results are 
shown by the black full curve, the D-RPA results by the red dot-dashed lines, P-RPA by the blue dashed lines, and the SC-RPA 
by the green dot-dashed lines. 



Fig. [5j for q = n/2, we see that the lowest state is quite well described by D-RPA as well as by P-RPA indicating 
again that quantum fluctuations are not very pertinent. On the other hand, the second excited state is not well 
reproduced by any one of the four approximations. Presumably this state has strong admixures of the multi-particle 
multi-hole type. In Fig. for q = 37r/4, we see for the first excited state a clear improvement of P-RPA over D-RPA. 
Again the other two states are not very well reproduced by any of the four approaches. 

Because of its particular structure, we split the q = it excited state case into two figures. We first compare, in 
Fig. [71 the S-RPA and the D-RPA results with the exact one. It is clear that there is no tendency for the D-RPA to 
reproduce the lowest exact level, which means that this level is not of the particle-hole type but of a different nature. 
Indeed, it is natural to assume that the existence of two HF states is responsable for a tunneling mechanism which 
splits the groundstate into a nearly degenerate doublet. In principle, above g c ~ 0.7, S-RPA looses its meaning and 
one might tend to ignore it. However, interestingly enough, the lowest real S-RPA level, which becomes complex for 
g > 1.8, reappears under the guise of the lowest level of the P-RPA, as Fig. [5J shows. We also show in Fig. [JJthc 
results of SCRPA which performs remarkably well in the spherical region (and beyond) for all four states. 

In Fig. [51 which also refers to q — w, we observe that the results of P-RPA improve the corresponding ones for 
the D-RPA for the two lowest states. However, for the two upper states, the situation is just the inverse: the P- 
RPA results are worse than the D-RPA ones. In this respect we should be aware of the well known fact that RPA 
overestimates the correlations and, therefore with attractive interaction the states come at lower energies than with 
an improved RPA. One sees this clearly for the two lowest states in Fig. 8 where the additional repulsion coming 
from the projection pushes the states upwards. The same upward push happens for the two upper states and the 
reasonable agreement of the D-RPA results with the exact solution must be interpreted as just an accident due to 
too much correlations in D-RPA. In the deformed region both D-RPA and P-RPA strongly deviate from the exact 
solution indicating that higher than ph-correlations are contained in those states in the deformed region. For g > 1.7, 
the lowest exact state is clearly not of the ph type. 

It is remarkable that the P-RPA leads to a level at almost zero excitation, for a surprisingly extended range of g 
values, mimicking a low lying state with odd parity which quickly becomes nearly degenerate with the groundstate. 
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Figure 7: (Color online) Eight site system. RPA frequencies vs g for the momentum transfer q — n. The exact results are 
shown by the black full curve, the S-RPA results by the red dotted lines, the D-RPA results by the red dot-dashed lines, the 
SCRPA results by the green dashed curves. 

In the next subsection it is shown that this level is very well reproduced by the Peierls-Yoccoz projection method in 
the version variation after projection @ . 

Let us finally mention that the states VpjV'gi * = 1, 2, used for constructing the RPA matrices J\f,T,V,J\f',T',V, 
were of the particle-hole type, that is, either pi G D, qj ^Dor pi ^ D, qj e D. 

C. The lowest excited state and the Peierls-Yoccoz method 

The Peierls-Yoccoz projection method has been devised to describe the rotational spectrum of a deformed nucleus. 
Let us assume that a prolate ellipsoidal nucleus rotates around an axis perpendicular to the symmetry axis. Let 
| denote the wave function of the rotated nucleus by an angle ft. According to the Peierls-Yoccoz method, the 
wave function of the nucleus is a superposition of the wave functions |$(f2)), 




dn/(n)|*(n)), 



(see Ref. [6]). From the symmetry of the nuclear Hamiltonian, it follows that /(O) = exp(imfi). For m = 0, a 
description of the groundstate is obtained. For m — 1,2,... a rotational band arises. In the Heisenberg model, the 
analogous symmetry to rotational symmetry is reflection symmetry, which is broken by the deformed HF state, so 
that we are dealing with a discrete group and the integral is replaced by a sum |W) = (/o + fiU\)\<&). From the 
symmetry properties of the Hamiltonian it follows that either f = fi (groundstate, with even parity) or / = — /i 
(excited state with odd parity). In this model, instead of a rotational spectrum, we find a doublet of states. We 
consider the variation after projection version of the method @. 




0,5 1 1,5 2 2,5 

g 



Figure 8: (Color online) RPA frequencies vs g for the momentum transfer q — ir. The exact results are shown by the full curve 
in black color. The performances of D-RPA, P-RPA and SC-RPA are compared. 

We have already discussed the ratio 

($|1 + E7i|#) 

It is, therefore, suggestive that we also consider the ratio 

which may describe the lowest state of odd parity. If we plot the difference 

5 = (gjg(l - E/QUl) + t/i)|$) 

(<l|l - C/i|$) ($|l + C/i|$) 

versus the coupling constant 5, the curve which is obtained falls almost on top of the lowest exact excited state in fig. 
6. In table I, we give S as a function of the coupling strength g, for the eight sites system. 

The feature that both states get degenerate in the infinite g limit is reproduced by the deformed HF state. Indeed, 
if we denote by <jy k the creation operator of a fermion in site k (k £ {1, 2, 3, 4, 5, 6, 7, 8}), then the groundstate wave 
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Table I: Description of the lowest excited state by the P-RPA and by the Peierls-Yoccoz projection 
ui (Exact) an (P-RPA) 8 



0.0 0.765367 0.765367 0.765367 

0.5 0.672780 0.636131 0.666271 

1.0 0.522674 0.568805 0.504778 

1.5 0.375722 0.460095 0.350829 

2.0 0.258764 0.000000 0.233469 

2.5 0.176111 0.000000 0.154336 

3.0 0.121093 0.000000 0.103791 



function is two-fold degenerate and reads <p\4>\<j)\(j)\\Q) or <f>\<f>\(j)\4>\\Q) . Now, in momentum representation we have 

4>\ = (Vie li ^ + ip-ie~ l! ^ + ip^e 1 ^ + ipse - * 3 **' 

= (e**(^ + i>- 7 e- ikn ) + e-**(V»-i + ^e 4fc ") 
+ c 1 ^ (^ 3 + 4>- 5 c- lk7r ) + c-^ (^_3 + ^ ik *)) , 

implying that the deformed HF becomes exact in the limit g — > oo. In this limit, one can neglect the kinetic energy 
in (1) and one may check that the above state is eigenstate to the interaction part of (1) alone. 



IV. CONCLUSIONS, DISCUSSION, AND OUTLOOK 



We have investigated the 8 sites and the 4 sites ID AFH model with HF and RPA type of theories. One may 
wonder why mean field type of theories should be applied to ID systems. However, the ID character only becomes 
specific getting close to the thermodynamic limit. As long as the dimensions of the discrete RPA problem stay 
rather low, i.e., the number of sites relatively small, whether one works in 3D or in ID has no real significance. This 
can, e.g., be seen in the case of other ID lattice models of the Hubbard type where RPA kind of approaches give 
very encouraging results, see for example [13, [H|. As is well known RPA may signal instabilities of the system. In 
order to avoid unstable RPA modes for large values of the anisotropy factor g, the symmetry unbroken HF state, 
i.e. the translational invariant HF state, defined in terms of a symmetrical single particle basis, should and has been 
replaced by a symmetry-broken HF state, defined in terms of a deformed single particle basis. Since the deformed 
HF state is degenerate, it is natural to apply the Peierls-Yoccoz projection method, in the version variation after 
projection Q. When the projection method is used, the groundstate is described as a superposition of two linearly 
independent deformed HF states and the deformed HF basis is shown to play a role even for small g values, for which 
the symmetrical HF state is stable. A formalism of the RPA type based on the projected groundstate, was outlined. 
The projected RPA has in common with the self-consistent RPA (SCRPA) [8j the fact that it is not related to a 
specific HF state, but rather to a correlated groundstate. However, the similarity between that theory and the present 
approach may stop here. We have found that the projection method considerably improves the HF treatment mostly 
as far as the groundstate energy is concerned, but also with respect to the lower RPA energies. While in D-RPA low 
lying levels are absent for g > 0.8, in P-RPA a low lying level which quickly approaches when g — ¥ oo, is found. 
The existence of this level is interesting, even though it deviates from the actual first excited state, which becomes 
degenerate with the ground state in g — > oo limit. Tunneling processes, which become inhibited in this limit, are 
poorly described by P-RPA and are ignored by D-RPA. We also observe that the projection method results become 
exact for the 4-site AFHM. In the 8-site problem, we applied the projection technique for both ground states, i.e. for 
positive parity and negative parity obtaining thus the first excited state with high accuracy. 

Concerning future developments, it might be interesting to formulate SCRPA also with a deformed single particle 
basis and apply the projection technique on top of it. Whether such an extension of the present approach will be 
feasible remains to be seen. 
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Appendix A 



As an example, we present some computational details for the 4-site system. 

The HF state reads |$) = ^V'ljO), hi = §7r, h 2 = -§7r. The deformed HF state 

\®)=(L(l\°)> 4 = - i = 1,2, a,/3eR, a 2 + f3 2 = 1, 

is determined by minimizing, with respect to /3, the expectation value 

($|ff|*) = -1 - {a 2 - [3 2 )V2 + |(a 2 - (3 2 ) 2 . 
The variation after projection procedure is implemented when we minimize, with respect to /3, the expression 

(gr#) -2(a 2 -(3 2 )V2 + g(a 2 -(3 2 ) 2 P 

£pr °^ wnm = ~ 1+ — i+{* 2 -p 2 ) 2 — ' Ui=e • (14) 

The minimum of £p ro j occurs for 



a 2 -(3 2 



V8 



For momentum transfer ±7r/2, we find 



A' Ur I 2 : (I { \) iVn 



T 



-V2V-/3 2 )(J J) (16) 

v = -1(1,-1) ^ 



For the projected state we still need the matrices Af',T',V. We find N' = AT, V = T, V = V. 
For momentum transfer 7T, we find 

H= {a 2 -/3 2 )diag(-l, -1,1,1) (18) 
r=V2(a 2 -/3 2 )diag(l, 1,1,1) (19) 
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V = 



/ 2a 2 /? 2 -l(a 2 -(3 2 ) 2 ~2a 2 /3 2 \{a 2 ~ fi 2 ) 2 \ 
I(a 2 -/3 2 ) 2 2a 2 /? 2 i( a 2 -/3 2 ) 2 _2a 2 /3 2 

-2a 2 /3 2 I(a 2 -/? 2 ) 2 2a 2 /? 2 -I( a 2 -,3 2 ) 2 

J 



(20) 



\i(a 2 -/? 2 ) 2 -2a 2 /3 2 -i(a 2 -/3 2 ) 2 2a 2 /3 2 
For the projected state RPA we still need the matrices TV', T', V. We find TV = TV, T = T' and 



/ 



V' = 



2 



_L 

2 



_ _L 

2 

Q 

_L 

2 









(21) 



■5 0/ 

The RPA frequencies based on the projected state are the roots of 

det [(TV + M')ui - (T + T') - g(V + V')] 
The parameter /3 should be such that it minimzes Eq given by (HU). 



Appendix-B 

We illustrate the determination of the matrices TV, T, V,TV, T', V for the 4 sites case. The fermion operators read 
V , azL>'0— 3»,"0s j but, for simplicity of notation, we write r/fcj, ^—3, ^1, ^>-i- The deformed fermion operators read 

£3 = aV>3 + M->-i, C-3 = a ^-3 + Ci — QT0i — /3ip—3, C-i = onj)-\ — /3"03- F° r momentum transfer q = 7r/2, the 
matrix TV reads 

\f — ( ^V,3;l,3 TVl^jS,! 

92P2 ; P191 - {$\0 N (q 2 p2]Piqi)\$), with (g 2 p 2 ), (p l9l ) e {(1, 3), (3, 1)}, being |$) = c£cl 3 |0>. The matrix TV 



where TV 
reads 



TV 



-^T,3;l,3 -^1,3;3,1 
■^3.1:1,3 -^3,1:3,1 



where TV^ 2P2;pi(?1 = (&\0 n (q 2 P2', Piqi)Ui\&) . The corresponding matrices T, V, T 7 , V" are similarly constructed. 
For momentum transfer q = n, the matrix TV reads 

TV = 



/M,-3;l 


-3 


M.-3; 


-1,3 


M.-3; 


-3,1 


M,-3;3, 




TV_1,3;1 


-3 




-1,3 


/V_ li3; 


-3,1 


N-l,3;3, 




TV_ 3 ,1;1 


-3 


/V_ 3 , 1; 


-1,3 


/V_ 3 ,l ; 


-3,1 


N-3,i;3, 




W3,-l ;1 


-3 


/v 3 ,_ 1; 


-1,3 


A/s.-i; 


-3,1 


A/3,-l;3, 


-1 



where TV 
reads 



<?2P2;pigi 



where TV' 



= ($|OAr(q 2 p 2 ;pigi)|$), with (q 2 p 2 ), ( Pm ) G {(1, -3), (-1, 3), (-3, 1), (3, -1)}. The matrix TV 



TV 



(^ICat (92^2 ; Pi9i )C^i I - The corresponding matrices T, V, T', V are similarly constructed. 





-3 


^i',-3; 


-1,3 


^{,-3; 


-3,1 


^-i,3;i, 


-3 




-1,3 


^1,3; 


-3,1 




-3 


•^-3,1; 


-1,3 


M-3,1; 


-3,1 




-3 


^3,-li 


-1,3 




-3,1 



Next, we exemplify the computation of the entries Vi t -3-i t -3 and V[ _ 3 . 1 _ 3 . From ([7]) we find 
Oy(l,-3;l,-3) = 4^3,1,-3,1^3^^-3 



- 4 V 3 ,-3,1,-1 ^3^13-0-1^1 - 4 W_3 -l,3,lV , i 3 V'ilV'l'03 
+ (2 ~ 2 W_ 3 _i,_3 ,-l) (V , - 3 V'-l'0-l'0-3 - V'lV'-lV'-l^l) 

+ (2 W-3,3,-3,3 - 2 ^1,3,1,3) (VwlV^l - ^13^3^3-0-3) ■ 



',t „/,t 



For instance, for the expectation value ($|?/;iiV ; -3V'3' ! /'i| < I ) ) = (OIC3C— 3 1 / 1 — l^LaV'sV'iC— 3C3 10) we obtain 

(0|C_3^ 3 |0)(0|C3^ 1 |0)(0|^ 1 d 3 |0)(0|^3C3 t |0) = ~a 2 f3 2 . 

On the other hand, the expectation value (§\ip j< _ 1 i> j L 3 i>3ipiUi\$) = (OIC3C— sV 7 — sV's'^iC— 3C^ |0) becomes 

<0|^l x |0> <0|C-3^I 8 |0> (Ol^xC'lalO) (Ol^C'slO) = a 2 /? 2 . 

Finally, ($|O y (l, 3; 1, 3)|$) = 2a 2 /3 2 , while 3; 1, 3)t/i|$) =0 



